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The very low temperature properties of two pyrochlore compounds, 
Yb 2 Ti 2 07 and Gd 2 Sn 2 07, were investigated using an ensemble of micro- 
scopic and bulk techniques. In both compounds, a first order transition is 
evidenced, as well as spin dynamics persisting down to the 20 mK range. 
The transition however has a quite different character in the two materials: 
whereas that in Gd 2 Sn 2 07 (at IK) is a magnetic transition towards long 
range order, that in Yb 2 Ti 2 07 (at 0.24 K) is reminiscent of the liquid-gas 
transition, in the sense that it involves a 4 orders of magnitude drop of the 
spin fluctuation frequency; furthermore, no long range order is observed. 
These unusual features we attribute to frustration of the antiferromagnetic 
exchange interaction in the pyrochlore lattice. 

PACS numbers: 75.40.-s, 75.25,+z, 76.75.+i, 76.80. +y 

1. Introduction 

Frustration of exchange interactions among magnetic ions can arise in 
different ways: the first to be studied was frustration induced by crystallo- 
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graphic disorder, leading to the so-called spin glass materials. Frustration 
can also appear in fully crystallographically ordered materials, if the geom- 
etry of the lattice is such that it prevents all pairs of exchange bonds from 
being satisfied throughout the lattice. The simplest example is the bidimen- 
sional triangular lattice with isotropic (Heisenberg) antiferromagnetic (AF) 
nearest neighbour exchange. Villain drew attention to geometrically frus- 
trated systems and showed that no Neel order can occur in a Heisenberg 
antiferromagnet on a three-dimensional lattice of corner sharing tetrahedra. 
The ground state of such a system was called a "cooperative paramag- 
net" , or a "spin liquid" state, where the spins undergo short range dynamic 
correlations down to T = 0. In the past decade, investigations of geometri- 
cally frustrated systems have developed to a large extent |2j , and a number 
of lattices have been shown to be prone to frustration: the bidimensional 
kagome lattice, made of corner-sharing coplanar triangles, the pyrochlore 
lattice, made of corner-sharing tetrahedra, and the garnet lattice, made of 
non-coplanar corner sharing triangles. 



Fig. 1. The pyrochlore lattice of R2M2O7 materials, where only the R sublattice 
is shown (black dots); for each i? 3+ ion, the local symmetry is threefold, with the 
appropriate local [111] axis as symmetry axis. 

The formation of the ground state and the low temperature spin dynam- 
ics have been particularly studied for the case of the Heisenberg antiferro- 
magnet in a pyrochlore lattice with nearest neighbour exchange [3]. The 
ground state configuration is shown to have a large degeneracy, stemming 
from the energy mimizing condition: J2i Si = 0, where the sum runs over 
the spins Sj of a "plaquette" or a "unit", i.e. either a triangle or a tetra- 
hedron. These ground states are not separated by energy barriers, and this 
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can have important implications as regards the spin dynamics: zero energy 
local or extended soft modes are possible, even at very low temperature. 
Another consequence is that any perturbation, like the dipole-dipole inter- 
action, next nearest neighbour exchange, ionic or exchange anisotropy, can 
select a particular ground state, resulting in a transition at a finite tem- 
perature. The perturbations can also alter the flat energy landscape and 
create energy barriers between ground configurations, which suggests that 
real frustrated pyrochlore systems can retain some features of spin-glasses. 

Compounds with formula R2M2O7, where R is a rare earth and M a 
transition or sp metal, crystallise in a structure where both R and M ions 
are located at the vertices of two interpenetrating corner sharing tetrahedra, 
or pyrochlore, lattices (FigQ). In the course of our study of rare earth py- 
rochlores, we have discovered novel behaviours in two materials, Yb2Ti207 
and Gd2Sn2C>7, which are the subject of this report. In these compounds, 
the rare earth alone is magnetic and it is located at a site with threefold 
symmetry (D^)- The paramagnetic Yb 3+ (4/ 13 ) and Gd 3+ (4/ 7 ) ions have 
distinct single ion properties: the crystal electric field splits the ground 
spin- orbit multiplet {J = 7/2} of the Yb 3+ ion into 4 Kramers doublets, 
whereas it has practically no influence on the Gd 3+ ion which has L = and 
S = 7/2. This implies that the ground Yb 3+ Kramers doublet, which alone 
is populated at low temperature, is described by an effective spin 1/2 with 
a g-tensor having uniaxial anisotropy, whereas the Gd 3+ ion is isotropic and 
has a g- factor very close to 2. 

Our investigations were carried out using local techniques: Mossbauer 
spectroscopy on the isotopes 170 Yb and 155 Gd, Muon Spin Relaxation (//SR) 
spectroscopy, Perturbed Angular Correlations (PAC) on the isotope 172 Yb 
and neutron diffraction on the one hand, and bulk measurements: magnetic 
susceptibility and specific heat on the other hand. Most of these measure- 
ments were performed down to the 20 mK temperature range. 

2. Yb 2 Ti 2 7 

In Yb2Ti207, a sharp peak in the specific heat had been evidenced near 
0.24 K 0] (Fig(2]left), but the nature of the transition had not been further 
investigated. A broad feature, peaking around 2-3 K, is also visible. As will 
be shown below, the ground crystal field doublet of the Yb 3+ ion is well 
isolated from the excited states. Therefore, the broad feature cannot be 
assigned to a crystal field Schottky anomaly, but it must be attributed to 
exchange spin correlations. The entropy released at the transition amounts 
to only 20% of the total magnetic entropy associated with a doublet, it! In 2. 
This behaviour shows that the correlations develop far above the tempera- 
ture of the transition (0.24 K), which is a characteristic feature of frustrated 
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systems |2j. 

In order to assess the crystal field splitting in this material, we performed 
PAC experiments on the isotope 172 Yb, which allow the thermal variation of 
the electric field gradient (EFG) tensor at the nucleus site to be measured. 
In axial symmetry, the EFG tensor is specified by its principal component 
Vzz = %^r( r = 0) alone, which is the sum of two terms: 

V ZZ (T) = V l z f + B Q (3J 2 Z - J(J + l)) r . (1) 

The first term is a temperature independent lattice charge contribution, and 
the second term is proportional to the 4/ shell quadrupole moment, which 
is temperature dependent due to the progressive population of the crystal 
field levels as temperature increases. 



50 




Fig. 2. In Yb^T^hO-j: Left: Specific heat (from Ref. 'S), and Right: thermal varia- 
tion of the principal component V zz of the electric held gradient at the nucleus site 
as measured by PAC on 172 Yb (solid symbols) and 170 Yb Mossbauer spectroscopy 
(open symbol); the line is a fit to a crystal electric Geld model (see text). 



The quantity V ZZ (T) decreases monotonously as temperature increases 
(FigHJright), and tends to V zz u at high temperature (~ —200 V/A 2 ) because 
the 4/ shell term vanishes when all the crystal field states are equipopu- 
lated. One observes that V zz (0) is quite small, which is to a coincidental 
cancellation of the lattice and saturated 4/ contributions, which are of op- 
posite signs and of comparable magnitudes. We also performed Mossbauer 
experiments at 4.2 K on 170 Yb diluted in non-magnetic Y2Ti20"7, which 
enabled us to measure the axially symmetric hyperfine tensor associated 
with the ground crystal field doublet. This tensor is proportional to the 
spectroscopic g-tensor, which is found to have components: g z ~ 1.80 and 
gj_ ~ 4.27. This shows that the plane perpendicular to the local [111] axis is 
an easy magnetic plane for the Yb 3+ ion at low temperature. From both the 
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V zz (T) thermal variation and the measurement of the ground state g-tensor 
(assuming that the latter is identical for Yb in Y2Ti207 and in Yt^l^Oy), 
we could determine the full Yb 3+ crystal field level scheme in Yl^T^Oy. 
In particular, the 3 excited doublets lie at 620, 740 and 940 K above the 
ground state. At low temperature, the physics is therefore governed by the 
anisotropic ground Kramers doublet alone. 




0.1 0.2 0.3 0.4 0.5 _ 2 — ^ J J 2 
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Fig. 3. InYb2Ti207: Left: magnetic susceptibility below 0.5 K showing the two FC 
and ZFC branches (inset: inverse susceptibility up to 20K), and Right: selected 
170 Yb Mossbauer absorption spectra between 0.036 K and 4.2 K; the lines are Gts 
as explained in the text. 

The magnetic susceptibility x{T) follows a Curie- Weiss law below about 
50 K, with a small positive paramagnetic Curie temperature p ~ 0.8 K 
(inset of FigEJleft). This can be indicative of a ferromagnetic net exchange 
interaction of magnitude ca. 1 K, but it does not preclude the presence of 
larger nearest and next nearest neighbour couplings with opposite signs. 
This latter assumption seems to be supported by the neutron diffraction 
data on a single crystal, to be described below. At 0.24 K, the temperature 
of the specific heat peak, an anomaly is seen in x(T) (FigO left). Below 
0.24 K irreversibilities between the Field Cooled (FC) and Zero Field Cooled 
(ZFC) branches appear, which are reminiscent of a spin-glass behaviour. 

Selected Mossbauer absorption spectra on the isotope 170 Yb {I g =0, 
I e =2, Eq = 84.3 keV), represented in Fig El right, reveal an apparently 
standard first order magnetic transition at 0.24 K. At 0.036 K, a five-line 
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magnetic hyperfine pattern is observed with a hyperfine field H^f — 115 T 
corresponding to a saturated Yb 3+ moment of ca. 1.15 fiB- As temperature 
increases, the spectrum remains unchanged up to 0.22 K; then a broad sin- 
gle line grows superimposed on the five-line spectrum, and at 0.26 K the 
single line alone is left. In the temperature range 0.036 K < T < 0.26 K, 
the hyperfine field associated with the five-line spectrum remains essentially 
constant. This fact and the coexistence, in a small temperature region, of 
the five-line pattern and of the single line, is characteristic of a first order 
transition. In principle, for a first order transition, one phase transforms 
into the other at the transition temperature by absorbing or releasing latent 
heat. In real systems, however, there is a narrow distribution of transition 
temperatures among the crystallites; this leads to the coexistence of the two 
phases in a small region around the mean transition temperature. 

For an anisotropic Kramers doublet, the saturated moment value de- 
pends on the orientation of the applied or exchange field and, except when 
the field is along a principal direction of the g-tensor, the moment and 
the field are not collinear. The relationship linking the modulus m of the 
saturated moment and the angle 6 it makes with the local z-axis is jSJ: 

cos 2 6 = Km ' , 2 
r z — 1 

where m± = \g^B and r = g±/g z - In Yb2Ti207 below 0.24K, from the 
measured m value and known r value, the angle between the moment and 
the local [111] axis is:9 ~ 45°. So the Yb moment does not lie in the easy 
plane perpendicular to [111]; the angle <p the exchange field makes with the 
local z-axis can also be obtained through the relation: tan 6 = r 2 tan ip. 
One gets: (p ~ 10°. So the exchange field is almost directed along the 
threefold symmetry axis. At the present time, the origin of these angles 
cannot be reasonably interpretated, but they may provide clues for a future 
investigation of the spin configuration in Yb2Ti2C>7. 

The characteristic time scale in these experiments is the hyperfine Lar- 
mor period for 170 Yb, which is about 10 -9 s. This value is in fact the center 
of the 170 Yb "Mossbauer window" where the spin fluctuation frequency 
can be measured from the lineshape [Hj. This has two consequences in the 
present case. First, it cannot be said whether the hyperfine field spectra 
observed below 0.24 K correspond to a static long range order (LRO) or to 
dynamic short range correlations with fluctuation frequencies smaller than 
the lower limit v\ ~ 10 8 s _1 of the "Mossbauer window". This problem will 
be dealt with below when describing the neutron diffraction and /iSR ex- 
periments. Second, the spectra above 0.24 K show a single line whose width 
decreases as temperature increases: this is the fingerprint of the "extreme 
narrowing" regime, where the spin fluctuation frequency is much larger than 
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the hyperfine coupling and increases as temperature increases. The fluctua- 
tion frequency can be measured as long as it is below the upper limit of the 
"Mossbauer window" for 170 Yb, i.e. v u ~ 5 x 10 10 s -1 . A physical hypoth- 
esis about the fluctuations is also needed: for paramagnetic fluctuations, a 
lineshape of the type developed in Ref.jS] is adequate, but for hyperfine field 
fluctuations (I.e. correlated moments), a stochastic lineshape as described 
in Ref. 7 must be used. In Yb2Ti207, the neutron diffraction measurements 
on a single crystal, to be described below, show that spin correlations are 
present well above 0.25 K. So we interpret the spectra in the temperature 
range 0.25 K < T < 0.9 K in terms of hyperfine field fluctuations, where Hhf 
jumps isotropically at random with a characteristic frequency v (7j. Then 
the dynamical broadening of the Mossbauer spectrum in the "extreme nar- 
rowing" regime is given by: 

AT R = ^, (3) 

where the hyperfine frequency u^t is fijHhf, being the magnetic moment 
of the excited 170 Yb nuclear state. The fluctuation frequencies obtained 
thereby are reported in Fig|S]and will be discussed below, together with the 
/xSR data. 

A first answer to the question about the presence of LRO in Yb2Ti207 
comes from neutron diffraction experiments on a powder sample (performed 
at the "Laboratoire Leon Brillouin", Saclay), which reveal there are no 
magnetic Bragg peaks below 0.24 K. Therefore there is no LRO in Yb2Ti207 
below the temperature of the specific heat peak. 

An insight into the dynamics of the moments is further provided by 
/iSR measurements (whose characteristic time is about 10~ 6 s), performed 
at ISIS (Rutherford Appleton Laboratory, Chilton, England) and at PSI 
(Villigen, Switzerland) L 8 j. For all temperatures down to 0.275 K (see Fig|l] 
a), the time decay of the muon depolarisation has an exponential form: 

aP z {t) = a z exp(-X z t) + a bg , (4) 

where a^g = 0.065 is a time independent background contribution arising 
from muons stopping in the Ag sample holder, and a z ~ 0.17. These data 
were recorded in a small longitudinal magnetic field B ex t = 2 mT in order to 
suppress the contribution of the nuclear moments to the depolarisation. The 
relaxation rate \ z is approximately constant (~ 0.1 MHz) as temperature 
decreases down to a few Kelvin, then it rises rapidly. An exponential decay 
for the depolarisation is usually observed when the fluctuation frequency v 
of the electronic moments is much larger than the dipolar (and hyperfine) 
coupling A between the muon spin and the electronic spin: v S> A ("ex- 
treme narrowing" limit). In this limit and in the paramagnetic phase (i.e. 
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Fig. 4. In YbiTiiO'j: (a) fiSR depolarisation spectra, measured at ISIS, on each 
side of the transition temperature (0.24 K) and (b) detail of the depolarisation at 
0.2 K at short times, measured at PSI; the lines are fits as explained in the text. 



when the electronic moments are uncorrelated) , the following relationship 
links \ z and v: 

A 2 

A 2 = 2-^, (5) 

where A p is the root mean square deviation of the distribution of dipo- 
lar couplings experienced by the muon spin in its interstitial stopping site. 
The rise of \ z below about 2 K is thus indicative of a slowing down of the 
electronic fluctuations as the spin correlations develop. Below the tempera- 
ture of the specific heat peak (0.24 K), the shape of the muon depolarisation 
changes drastically (see the 0.2 K data in Fig01): it is no longer an exponen- 
tial function of time, but it shows a rapid depolarisation within about 0.2 [is 
followed by a slow quasi-exponential decay. This shape remains unchanged 
down to the lowest temperature of the experiment, 0.04 K. The first remark- 
able thing to notice is that no oscillatory signal is observed, although in this 
temperature range a hyperfine field is observed in the Mossbauer spectra. 
This means that no LRO is present, in agreement with the powder neu- 
tron diffraction data. The shape of the low temperature depolarisation can 
actually be accounted for by a dynamic Kubo-Toyabe decay (dashed line 
in Fig0J), except for a small discrepancy around 0.03 /is: the dip present 
in the theoretical curve is blurred in the experimental data. A dynamic 
Kubo-Toyabe lineshape describes the general case of a muon spin coupled 
to an ensemble of electronic spins with fluctuation frequency u, and with a 
root mean square isotropic coupling A. Its high frequency limit is an ex- 
ponential function with a relaxation rate given by expression ©. In order 
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to reproduce correctly the low temperature depolarisation, a solution con- 
sists in introducing a gaussian distribution of A values, which could be a 
means of taking into account the spin correlations 9 j. A good fit is thereby 
obtained (solid line in Fig@] b), with a rather low fluctuation frequency: 
v ~ 10 6 s -1 . The mean value of the muon - 4f spin coupling, in units of 
field, is: Alt/i^ — 5.7 mT, where 7^ = 851.6 Mrd s _1 T _1 is the muon gyro- 
magnetic ratio. Thus the presence, below 0.24 K, of a /^SR signal consisting 
of a dynamic Kubo-Toyabe decay and of a hyperfine field in the Mossbauer 
spectra, means that short range dynamic spin correlations are present down 
to 0.04 K. The measured fiSR fluctuation frequency of 10 6 s _1 is accordingly 
lower than the Mossbauer threshold v\ ~ 10 8 s _1 . 

Above the temperature of the specific heat peak, the electronic fluctua- 
tion frequencies v can be deduced from the X z values through formula (jSJ 
if the high temperature A p value is known and also, in principle, if the 
spin correlations are weak. Usually, in the paramagnetic phase, fluctuating 
moments of \\ib (the size of the moments in Yb2Ti2C>7) correspond to a 
mean dipolar field of the order of 50 mT. So the mean Alt value of 5.7 mT 
obtained below 0.24 K, where the spin correlations are well developed, is 
too small and cannot be used at high temperature where correlations are 
weaker. So in order to extract v from the A 2 data above 0.24 K, we adjusted 
A p in expression (jSJ) so that the frequency values match those measured by 
Mossbauer spectroscopy (see FigEJ). We obtained Ajjt — 80 mT, which is 
the correct order of magnitude. In the paramagnetic phase, Aht can be cal- 
culated for a known crystal structure by assuming that the muon stops in a 
given interstitial site (and assuming also a dipolar only interaction with the 
rare earth spin) ^Uj. In the present case of an axially symmetric Kramers 
doublet with g-tensor {g z ,g±}, Aht is given by the expression: 

*%T = ~ + fliEa], (6) 

8irgj 

where Si and E2 are lattice sums whose expressions are given in Ref. jlOj . 
There are three possible interstitial sites for the muon in the pyrochlore 
lattice for which: Ajjt/1u — 150 mT. This value is larger by a factor of 2 
than that obtained by scaling the ^SR data to the Mossbauer data. This is 
an acceptable agreement, if one keeps in mind that expression is valid 
for uncorrelated moments only. In fact, at least below 2 K, the specific heat 
data and the neutron diffraction data on a single crystal (see below) show 
that spin correlations are present. Therefore, the muon spin is likely to 
interact with a well defined rare earth moment, of modulus 1.15 [is (the low 
temperature value), within correlated clusters, rather than with a param- 
agnetic Kramers doublet. With this assumption, the standard calculation 
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of Aht due to randomly oriented moments yields Ajjt Hp, — 80 mT for the 
three sites, which is exactly the value derived above. 

The overall thermal variation of the Yb spin fluctuation frequency is 
represented in Fig|5] and constitutes the main result of our study in this 
compound: the fluctuation frequency undergoes a first order abrupt drop 
at the temperature of the specific heat anomaly (0.24 K), falling from the 
10 4 -10 5 MHz range to 1 MHz. So the transition involves the time domain and 
is reminiscent of the liquid-gas transition, which is first order and involves 
an abrupt drop of the mean collision frequency between atoms as one enters 
the liquid phase. 
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Fig. 5. Thermal variation of the Yb spin fluctuation frequency as measured by 
/iSR (solid circles) and 170 Yb Mossbauer spectroscopy (open squares); the dashed 
line is the lower limit vi of the Mossbauer window for measurement of fluctuation 
frequencies. 

The low temperature phase in Yb2Ti20"7 is a "spin-liquid" phase, as 
predicted by Villain [Q for spins in a lattice of corner sharing tetrahedra, 
where spin dynamics is present down to the lowest temperature. However, 
this state is not reached "smoothly" by a continuous decrease of the fluctu- 
ation frequency, but by a first order transition in the fluctuation frequency. 
At this stage, it is not known whether the transition is accompanied by 
some change in the spatial spin correlations. The answer to this question 
comes from a neutron diffraction study of a single crystal, to be described 
next. 

The neutron diffraction experiments on a single crystal of Yb2Ti2C>7 
were performed at the Institut Laue Langevin (Grenoble, France), in the 
temperature range 0.04 K - 30 K. Like the previous neutron diffraction mea- 
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surements on a powder sample, no magnetic Bragg peaks were discovered as 
temperature decreases from 4.2 K to 0.04 K. But mapping of the magnetic 
elastic scattering in (110) planes reveals the presence of "diffraction rods" 
along the [111] direction in reciprocal space (see FigEJ). This is indicative of 
the presence of bidimensional antiferromagnetic spatial correlations, prob- 
ably within planes perpendicular to [111]. No anomaly is found at 0.24 K: 
the intensities of the peaks obtained from Q-scans across the rod (see Fig|7] 
left for the scan along line 1) increase steadily on cooling from about 25 K 
down to 0.04 K (see Fig|7]right, where the data below 1.4 K are not shown). 
The width of the peaks in Q-space does not vary with temperature and 
corresponds to a correlation length £ ~ 4 nm (about 4 times the parameter 
of the unit cell, which contains 16 Yb ions). The short range correlations 
therefore involve a few hundred Yb ions. 




0.5 1 1.5 2 2.5 

(0,k,k) 



Fig. 6. [color] Map of the magnetic neutron scattering in a single crystal of 
Yb2Ti20? in the (110) plane of the reciprocal space at 1.4K, showing the "diffrac- 
tion rod" along [111]. The red lines labelled 1 and 2 are the directions along which 
Q-scans were performed. 

The magnetic scattering intensity near Q=(1.5,1.5,1.5) starts to develop 
near 25 K, whereas that near Q=(1.9,1.9,1.9) starts to grow at a lower tem- 
perature (about 13 K). Both peaks however reach half their maximum in- 
tensity around 2-3 K. This suggests that these spin correlations are driven 
by the exchange evidenced in the specific heat bump, which has the same 
energy scale of 2-3 K. The short range bidimensional spin correlations in 
Yb2Ti2C>7 thus build up monotonically as temperature decreases; the range 
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Fig. 7. [color] Left: Q-scan across line 1 of the map in Fig\$at 1.4 and 27 K, and 
Right: thermal variations of the amplitudes of the peaks near (1.5,1.5,1.5) across 
line 1 and near (1.9,1.9,1.9) across line 2 of the map in FigQ3 



of these correlations, about 4nm, is temperature independent and in par- 
ticular shows no anomaly at 0.24 K. This reinforces the picture sketched 
above, of a transition involving the frequency domain only. 

3. Gd 2 Sn 2 7 

As shown in the inset of FigJSJleft, the inverse susceptibility in Gd 2 Sn 2 0y 
follows a Curie- Weiss law with a paramagnetic Curie temperature 6 P ~ 
— 10 K, indicative of an antiferromagnetic exchange interaction. As the 
Gd 3+ ion is isotropic, this compound can be therefore expected to be a 
good realisation of an AF Heisenberg frustrated system, with no Neel order 
down to T = 0. However, both the magnetic susceptibility and the specific 
heat, shown in Fig|S] right and left respectively, evidence an anomaly at 
1 K. The susceptibility x(T) presents furthermore a sizeable irreversibility 
between the FC and ZFC branches below 1 K. The anomaly of the spe- 
cific heat C P (T) at IK reaches the very large value of 120JK mol.Gd , 
whereas the expected jump at TV for a second order magnetic transition is: 
AC P = |-R s( |+ 1 )+o. 5 HU, i.e. 20.4JK- 1 mol.Gd- 1 for Gd 3+ with S=7/2. 

This high AC P value is indicative of a first order transition (the specific 
heat diverges in principle in this case), which will be confirmed by the 155 Gd 
Mossbauer measurements to be described below. Below 1 K, the specific 
heat varies as T 2 . Classically, the specific heat due to the thermal excitation 
of magnons in a 3-dimensional Heisenberg AF varies as T 3 , and as T 2 in 
the bidimensional case. So the observed quadratic variation in Gd 2 Sn 2 07 
suggests a bidimensional AF structure. The inset of FigJSl left shows the 
magnetic entropy variation as temperature increases; at the transition, only 
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Temperature (K) Temperature (K) 

Fig. 8. In GC/2S112O7: Left: magnetic susceptibility measured with lmT, with the 
Field Cooled (FC) and Zero Field Cooled (ZFC) branches; inset: inverse suscep- 
tibility below 30 K, and Right: thermal variation of the specihc heat; the dashed 
line is a T 2 law; inset: thermal variation of the entropy. 



40% of the total entropy Rln8 = 17.3 JK _1 mol.Gd _1 has been released, 
and the full paramagnetic degrees of freedom are recovered only at 8-10 K. 
As in the case of Yb2Ti207, this is due to the presence of short range 
order developing well above the transition temperature and it evidences the 
presence of frustration. 

Selected Mossbauer absorption spectra on the isotope 155 Gd (7 s =3/2, 
I e =5/2, £0=86.5 keV) in Gd2Sn20"7 are shown in Fig|5J A clear change in 
the lineshape is observed between 1.05 and 1.1 K (see Fig|H]left): at 1.1 K 
and above, the spectrum is a pure quadrupolar hyperfine pattern character- 
istic of the paramagnetic phase, while at 1.05 K a magnetic hyperfine field 
has appeared, evidencing short or long range magnetic order. As the tem- 
perature is further decreased down to 0.027 K, the hyperfine field increases 
slightly to reach a saturated value of 30 T. This behaviour is the hallmark 
of a first order transition, as also inferred from the height of the specific 
heat peak; in Gd2Sn207 however, the coexistence region between the two 
phases could not be detected like in Yb2Ti20y. The fits furthermore show 
that each Gd 3+ moment is perpendicular to the [111] local symmetry axis. 

The question arises here also as to whether there is LRO below 1 K, 
or whether the correlated moments fluctuate with a frequency lower than 
v\ ~ 3 x 10 7 s _1 , which is the lower limit of the "Mossbauer window" for 
155 Gd. The answer will be provided by the fiSK measurements, but before 
we turn to their description, we will examine the problem posed by the 
lowest temperature Mossbauer spectrum. 

The spectrum at 0.027 K (see Fig E] right) indeed deserves particular at- 
tention, because the temperature which can be measured from the spectrum 



14 



ZAKO03 PRINTED ON FEBRUARY 2, 2008 




-6 



-3 



) 
v (mm/s) 



3 



6 



-6-3 3 



6 



velocity (mm/s) 



Fig. 9. 155 Gd Mossbauer spectra in Gd2Sii20j: Left: at selected temperatures 
between 0.027K and 4.2K, and Right: at 0.027K; the lines are theoretical curves 
for various values of the hyperfine temperature T^f. dashed line, T^j —0.027 K, 
thick solid line, which fits the data well, T^t —0.090 K, and thin solid line: high 
temperature limit, Thf > 0.2 K. 



(~ 0.090 K, referred to hereafter as the hyperfine temperature Thf), is much 
higher than that of the sample (given by the thermometer): ~ 0.03 K. The 
possibility of measuring the absolute temperature from the relative line in- 
tensities of the Mossbauer spectrum has been since long recognised |12j . The 
conditions are that the spin I g of the ground nuclear state is non-zero, and 
that the hyperfine splittings are of the same order of magnitude as k^T . 
Then each resonant transition has an intensity which is weighed by the 
population of the ground hyperfine level it arises from. The temperature of 
the hyperfine levels (in principle equal to the lattice temperature) can then 
be obtained directly from the spectrum. This procedure is impossible for 
170 Yb, which has I g =0, but should work for 155 Gd, which has I g =3/2. In 
Gd2Sn207, the combined quadrupolar and magnetic hyperfine interactions 
yield four ground hyperfine levels with energies 0, 0.05, 12.1 and 15.9 mK, 
i.e. two quasi-degenerate doublets separated by a mean hyperfine split- 
ting Ahf — 14 mK. Then the hyperfine temperature can be measured if it is 
below about 0.15 K. On Fig^Jright are shown, together with the experimen- 
tal spectrum at 27 mK, the theoretical spectra expected for Thf = 27 mK 
(dashed line), for Thf > 0.2 K (thin solid line) and for Thf = 90 mK (thick 
solid line). This last line goes perfectly well through the data points, and 
it is then clear that the hyperfine levels are "hotter" than the lattice, i.e. 
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they are out of thermal equilibrium. 

0.65 i . 




Time (a.u.) 

Fig. 10. Temporal evolution of the ground state population pe(t) of a (nuclear) spin 
1/2 doublet, as the (hyperhne) Geld reverses randomly in time, at arbitrarily chosen 
instants 5, 9, 17 and 20. The ratio of the doublet splitting to ksT is chosen to 
be about 2, so that the Boltzmann population of the ground state is 0.63. Dashed 
line: T\ <C r: thermalisation is effective; dash-dotted line: T\ 3> r: thcrmalisation 
is impossible, and solid line: T\ ~ r. 



The interpretation we give for this phenomenon is the following |13| : 
the hyperfine levels have no time to reach thermal equilibrium, with a time 
constant T\, because the hyperfine field reverses with a time constant r 
shorter than or of the same magnitude as T\. This is illustrated in Fie llOl 
which depicts the temporal evolution of the population pi of the ground 
state of a fictitious spin 1/2 (modelling the hyperfine levels) as a function 
of the ratio T\jr. The steady state population (i.e. time averaged) (pe) 
is the Boltzmann population pf if hyperfine relaxation occurs very rapidly 
with respect to r (dashed line in FigllUj). it is 0.5 (equipopulation) in the 
reverse case (dash-dotted line in Fig llOj) . and it should be a function of the 
ratio fi = Ti/t if the latter is of the order unity (solid line in Fig llOj) . In 
order to obtain this function, we devised a model of a spin 1/2 submitted 
to a magnetic field of constant magnitude, but which undergoes "flips" at 
random instants of time. The model provided an analytical solution for the 
steady state population of the ground level: 



(Pi) 



1, 1 

- 1 + 

2 y l + 2/i 



tanh • 



A , 
2k B T' 
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16 



ZAKO03 PRINTED ON FEBRUARY 2, 2008 



where A is the Zeeman splitting of the spin 1/2 doublet. Assigning to 
the doublet an effective (in our case hyperfine) temperature (this is always 
possible for a spin 1/2 system), we obtain: T^f — T(l + 2/x) in the range 
JzbT > 2A which is the temperature range of our experiments. Applying 
this relation to the case of Gd2Sn2C>7 yields \x ~ 0.8. 

The fact that the hyperfine populations are out of equilibrium at 0.027 K 
reveals therefore the presence of two forms of dynamics of the Gd 3+ elec- 
tronic spin. The first one consists in spin flips between one (or more) direc- 
tion^) which, as the Gd 3+ moment is proportional to the hyperfine field, 
correspond to flips of the hyperfine field sensed by the 155 Gd nuclear spin. 
Besides, the nuclear (hyperfine) relaxation, at such low temperature, can- 
not be due to phonons. The only plausible mechanism is nuclear relaxation 
driven by coupling to electronic spin- waves ^H, which are thus the second 
form of spin dynamics evidenced at 0.027 K in Gd2Sn207. The presence of 
spin-waves at very low temperature is also demonstrated by the /aSR exper- 
iments, to be described below. The probe in that case is not the nuclear 
spin of the rare earth atom itself, but the spin of positive muon lying at 
an interstitial site; the relaxation mechanisms are however similar in both 
cases. 

The /tSR experiments were performed in zero magnetic field between 
20 mK and 100 K at PSI. Above 1 K, the decay has an exponential form, 
and the relaxation rate \ z has a temperature independent value of 2 MHz. 
Below 1 K, the decay signal changes and presents an oscillating component, 
indicative of magnetic LRO. The oscillations are clearly visible only at very 
short times (see Fig llll left). 



0.05 




Time (/as) Time (us) 

Fig. 11. Time variation of the fiSR depolarisation in Gt^S-i^O^ at 20 mK, measured 
at PSI: Left: short time part showing the rapidly oscillating components, and 
Right: long time part with a 100 channel binning in order to integrate out the 
oscillations. 
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An exponential damping component is also present, which is best visible 
when the counting channels are grouped into batches of 100 ("binning" of 
100 channels), as shown in Fig llll right where the line is a fit to an exponen- 
tial function. Usually, in the magnetically ordered phase, the depolarisation 
in a polycrystalline sample is well described by the following expression (as- 
suming a single stopping site for the muon) [Tp] : 

aP z (t) = — [exp(— X z t) + 2 exp(— Aj_i) cos u^t]. (8) 
3 

This is equivalent to a depolarisation function where one third of the im- 
pinging muons have their spin parallel to the internal field (first term) and 
two thirds have their spin precessing in the internal (dipolar) field (second 
term). In expression a is the amplitude (or asymmetry) of the depolar- 
isation, which in principle does not change when the temperature is lowered 
through the transition; X z is the longitudinal relaxation rate, due to the 
coupling of the muon spin with the electronic degrees of freedom (in the 
LRO phase, the spin- waves) ^H]; oj^ = 'y^Bi is the Larmor frequency of the 
muon spin in the internal field Bj, arising from the spontaneous electronic 
moment, and A^ is a damping factor which is generally taken to be due 
both to a distribution in £>j values and to fluctuations of Bj. As discussed 
above, the depolarisation in Gd2Sn207 below 1 K is described by expression 
(jSJ) in a first approximation. However, some features are not fully under- 
stood at the present time: the asymmetry a below 1 K is half the value 
above IK (a ~ 0.2); there are two oscillating components of equal weights, 
with dipolar fields of 206 and 441 mT, each with a sizeable phase shift, 
respectively 24 and 69°. A particular magnetic structure could be at the 
origin of these peculiarities. We will only mention here that the two dipolar 
fields are temperature independent between 20 mK and 0.7 K, the highest 
temperature where the oscillations are clearly visible. The dipolar fields at 
the muon sites are due to the spontaneous Gd 3+ moments, and their lack 
of thermal variation is in line with the Mossbauer spectroscopy finding of 
an essentially temperature independent moment value, i.e. of a first order 
transition. Their order of magnitude (200-400 mT) is in agreement with the 
rough rule stating that a moment of l/i_e gives rise to a dipolar field of 
50 mT, as the Gd 3+ saturated moment is 7 \ib- 

The thermal variation of the relaxation rate X z is shown in Fig |121 Be- 
low 1 K, it reflects the muon spin lattice relaxation through the coupling 
with spin- waves, and it is defined by Eqn.(JHJ); above IK, it reflects the fast 
relaxation due to coupling to the short range correlated, and eventually 
paramagnetic, Gd 3+ spins. 

The transition at 1 K is marked by a sharp anomaly in X z (T) . The 
remarkable feature is that \ z below 1 K shows no thermal dependence, re- 
maining at the constant value 0.6 MHz down to 20 mK. This is anomalous, 
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Fig. 12. Thermal variation of the muon longitudinal relaxation rate \ z in 
Gd 2 Sn 2 0j. Below IK, X z is defined by Eqn.^j; above IK, it is the relaxation 
rate of the observed exponential decay (the data for T > 4K were obtained with a 
longitudinal magnetic field of 10 mT). The dashed line at 0.6MHz underlines the 
constant relaxation rate below 1 K. 



as relaxation of the muon spin by spin-waves should vanish at T = 0. For 
example, it was shown that, for a Heisenberg ferromagnet, X z varies as T 2 to 
a good approximation if one considers two-magnon Raman processes, which 
are the most likely to occur for relaxation between the quasi-degenerate spin 
states of the muon [TS|. An analogous calculation for a Heisenberg bidimen- 
sional AF, which could apply to Gd2Sn207 as suggested by the specific heat 
data, yields a T 3 variation. None of these thermal laws matches with ex- 
periment. The explanation for the temperature independence of A 2 must 
probably be looked for in the unusual spin-wave spectrum of Gd2Sn207, 
where soft spin-wave modes with a finite density of states at zero energy 
could exist. A saturation of \ z below about 2K has also been observed 
in another pyrochlore material, Tb2Ti207 ^7j- But this compound shows 
no phase transition to an LRO phase at low temperature and remains a 
"correlated paramagnet" down to T = 0. 



4. Conclusions 

Both studied pyrochlore materials, Yb2Ti207 and Gd2Sn2C>7, present a 
low temperature state where spin fluctuations persist as T — > 0. 

In Yb2Ti20"7, this state seems to be a spin-liquid phase where short 
range dynamic spin correlations are present. Such a phase is expected for 
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geometrically frustrated isotropic systems with antiferromagnetic interac- 
tions. However, it is not reached by a smooth decrease of the fluctuation 
frequency as the temperature decreases, but through a first order transition 
(at 0.24 K) where the spin fluctuation frequency drops from the GHz to the 
MHz range. Below 0.24 K, this frequency could be measured directly by 
the ^SR experiments and it remains constant down to 40 mK. No anomaly 
could be detected at the transition for the spatial spin correlations: they 
build up from 10-20 K and remain short range, with a correlation length of 
4 nm, down to the lowest temperature. 

In Gd2Sn2C>7, the low temperature phase has long range magnetic order 
and it is reached through a standard first order transition. However, spin 
dynamics is present down to 30 mK: it was detected indirectly by 155 Gd 
Mossbauer spectroscopy through the finding of nuclear levels which are out 
of thermal equilibrium. This requires the presence of both electronic spin 
flips and of nuclear relaxation by coupling to spin waves, with similar char- 
acteristic time scales. The presence of spin waves down to 20 mK was also 
evidenced in the /iSR measurements by the non-vanishing relaxation rate of 
the fi + spin. 

These findings illustrate the unconventional dynamics occurring in the 
frustrated pyrochlores at very low temperature and the diversity of routes 
leading to a ground state where spin fluctuations persist as T — > 0. 
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